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Abstract 

We give a method for computing the iterated Laplace transform of the sojourn time in an union 
of intervals for hnear diffusion processes. This random variable comes from a model occurring in 
biology concerning the clustering of membrane receptors. The way used hinges on solving differential 
equations. We finally have a look on the particular case of Brownian motion and we provide a 
representation for the Laplace transform of its local time in a finite set. 
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1 Introduction 

We consider a diffusion process {Xt)t>o evolving on the real line M with infinite lifetime (that is without 
any killing). Let uj, wi, . . . , m„, w„ be real numbers such that wi < wi < • • • < u„ < i;„. Set also vq = — oo 
and Un+i — +00. Let us now introduce the union of intervals 

n 

together with the sojourn time that the process {Xt)t>o spends in up to a fixed time t: 

rt 

Of course, introducing the sojourn times of {Xt)t>o in [ui, u^], i £ {1, . . . , n}, 

Jo 

it is plain that time Tt can be decomposed into 

n 
i=l 

The functional Tt is one of the most classical and historical random variables introduced in stochastic 
processes theory. It has been of great interest for many researchers in particular when the set E is an 
infinite interval [u, +00) or (— 00, w]-cases leading to the famous arcsine law when {Xt)t>o is Brownian 
motion for instance-or a bounded interval [ujw], or when E is reduced to a point {u}. In this last case, 
it is related to the famous local time after some normalization. 
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The case of an union of intervals for E has appeared to us in a biological context that we describe here. 
Consider a cellular medium. The plasmic membrane is a place where many interactions occur between 
the cell and its direct external medium. In classical models, the membrane is described as a fluid mosaic 
made of different constituents; certain, the so-called ligands which possibly come from the extracellular 
environment, randomly roam along the membrane while other group into still clusters (the so-called rafts 
viewed as receptors) which form inhomogeneities on the membrane. Few functional properties of these 
inhomogeneities are known. A particular phenomenon is the binding mechanism between ligands and 
receptors. As a matter of fact, a quantity of great interest for biologists is the docking-time^ this is the 
time that ligands and receptors bind up to a finite time (the duration of the experience). Indeed, this 
quantity is an important indicator of affinity /sensitivity of the ligands for receptors. For a more accurate 
description of this biological context, we refer the reader to the paper by Care & Soula [2]- 

In [3] , we proposed a random walk model for simulating this problem: the membrane is viewed as the 
integer line Z (or a finite torus X/NZ,), the receptors are fixed numbers ai, . . . , in Z and the ligands 
evolve along the membrane like independent random walks (S'i(z))igN, (S'2(i))igN, ■ • • on Z. Hence, the 
movement of the family of £ ligands is modeled by the .^-dimensional random walk (S(i))igN defined, 
for any i S N, by S(i) = {Si{i), . . . , Si{i)). The docking-time is proportional to the sojourn time spent 
by the random walk (S(i))igN in the set uj^iZ-'"^ x {oi, . . . ,0^} x Z^~^. We provided a methodology 
for computing the probability distribution of this sojourn time. At the end of 3.^, we addressed the 
continuous counterpart to this problem: replace the ^-dimensional random walk by an ^-dimensional 
Brownian motion (or a more general diffusion process) and compute the distribution of the time spent 
by Brownian motion in a set of the form Uj^iW~'^ x E x M^^-' where i? is a set modeling rafts, that is, 
E is an union of intervals. 

In this paper, we consider the case of the dimension £ ~ 1; this is the evolution of one ligand roaming 
as a diffusion process along a linear membrane (viewed as R). In this case, the docking-time is quantified 
by the sojourn time of the diffusion process in an union of intervals, namely Tt. Our aim is to describe 
the probability distribution of Tt as well as the joint distribution of (Tt,Xf). For this, we provide a 
representation of the iterated Laplace transforms of these distributions which could be easily implemented 
on a computer. Actually, our method allows us to compute the joint distributions of {T^, . . . ,T") and 
{T^, . . . ,Tp,Xt). We finally have a look on the particular case of Brownian motion and we provide a 
representation for the Laplace transform of its local time in a finite set. 



2 Settings 

Let us introduce the transition densities of the process (^4)00: for any t > and x, y G 

p{t- X, y) = ¥,{Xt G dy}/dy - P{Xt G dj/j^o - x}/dy 

together with its A-potential: for any A > and x, y G M, 

1*00 

Px{x,y) ^ I er^^ p{t]x,y)At. 



Jo 

Let T> be the infinitesimal generator of the diffusion {Xt)t>o- This is a second-order differential 
operator. Since the lifetime of {Xt)t>o is infinite, the killing rate vanishes and then T> has the form 

The notation means that the operator V acts on the variable x. It is well-known that the density p 
solves the Kolniogorov backward equation 

dp 

V^p{t;x,y) = —{t;x,y), p(0+;x,y) = S{x - y) 
and that the potential solves 

1^xP^{x, y) = Xp^{x, y) - 5{x - y). 

The potential is the Green function of the operator I? — A id. 

Let us introduce a basis of fundamental solutions {c, d\ of the linear second-order differential equation 
2?u = Am such that c is increasing and d is decreasing. Then the potential admits the following 
representation (see, e.g., [Tj p. 19): 

Px{x, v) = -i-T c(min(a;, y)) d(max(x, y)) 
w{y) 



2 



where w{y) = [c'{y)d{y) - c{y)d' (y)]/ K{y) with K{y) = 2/a{y)'^. 

We see that the potential fulfihs the foUowing properties. The function x p^{x,y) is continuous 
on R and, in particular, at y. Moreover, it is diffcrentiable on ]R\{2/} and its derivative admits a jump at 
y. More precisely, we have the conditions below at y: 

8n dn (2-1) 

^{v^.v)-^{v-,y)^-<y)- 

We also need to introduce the first hitting time of any level u: 

Tu = inf{< >0:Xt^u}. 
The Laplace transform of the distribution of t„ is given by (1 p. 18) 

—— if x < u, 
(diu) 

In this paper we compute the iterated Laplace transform of the random vector — {T^^ . . . , T"): 

<^,»= / e-^*E4e-<'^^^'>)di 
Jo 

as well as that of the random vector (Tj, Xt): 

^.S^^y) = / e-^* [¥.,{e-<^'-'^^>,Xt e dy) /dy] dt. 

JO 

In the foregoing quantities, /ii is a vectorial argument: /x = (^i,...,^„) and < /x, Tf > is the inner 
product of the vectors /i. and T^: 

n 
i=l 

Notice that we use bold symbols for vectorial variables while unbold symbols are related to scalar variables. 
The iterated Laplace transform of the random variable Tf can be immediately deduced from ip^ {x) by 
identifying all the components of ^t: if /ii = • • • = /i„ = /i, then < ^, Tt >= /iTj. 

For this, we begin by writing integral equations satisfied by the functions ip^ ^ and ^ which lead to 
ordinary linear differential equations associated with some regularity conditions. Of course, we have 

/+00 
^a,m(^' y) dy- 

So, because of this relationship, it should be enough to compute ip^^^ and next integrate it with respect 
to y for deriving ip^ . Actually, this way appears to be untractable from the obtained formulas, so we 
shall evaluate separately both functions (p^ ^ and V'^ ■ 



3 Some equations 

In this section, we derive some integral equations (Subsection 13. ip and next deduce some differential 
equations (Subsection 13. 2|) satisfied by the functions ip^^^ and ^. The method we use is standard, this 
is the famous Feynman-Kac approach. Nonetheless, in order to facilitate the reading of the paper, we 
provide some details. Our results come from solving differential equations p.6p . p.7p below associated 
with the boundary values p.Sp . (|3.9p . p.lOp . The aforementioned equations as well as the additional 
conditions directly come from integral equations p.ip and p.2p . 
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3.1 Integral equations 

Proposition 3.1 The functions ip^ ^ and tp^ ^ solve the following integral equations: 



2^ ^ r'^i 

f^J^) = J^^f^^ PAx,z)ip^Jz)dz, (3.1) 



1=1 



AJx, y) = p^{x, y)-^fi, / p^{x, z)i^^Jz, y) dz. (3.2) 



Proof 

Let us introduce the quantity 

xSt;x,y)^E,{e~<''''^*>,Xt e dy)/dy. 

We write that 



1-e 



where {9s)s>o is the usual shift operator defined by XfO 9s = Xs+t for any s,t > 0. Next, we apply the 
Markov property: 

xSt;x,y) - pit;x,y) ^ ^E,{1 - c-<''''^^> ,Xt € dy) /dy 

= -Y.^'^ [E,(]l[„,.,,](X,)e-<'^'T'-=°^=>,Xt_,o0, edy)/dy]ds 

= ~J2^'^ / ^ "^^'^ [E,{c-<''''^'-^>,Xt-s e dy)/dy] ds 

i=l "'0 Jui 

= / P(«;^'^)x^,(^-s;^,y)dsdz. (3.3) 

By applying the Laplace transform to p.3p . we get Eq. p.2p which in turn yields Eq. p.ip by integrating 
with respect to y on R and using 



/+00 i>oo 
p,{x,y)dy^ / e-^'P4Xt G 



We also state the following result which will be useful. 
Proposition 3.2 For x £ (— oo,ui), we have 

cp.Jx) = IE4e-^^-)(^,,>i) -\)+\, (3.4) 
'^>.Jx, y) - p^{x, y) = E^ (e"^^"i ) [V';,,^(ui , y) - p^wi, y)] • (3.5) 

Proof 

Suppose X < ui. When t^-^ > t, the process (Xs)s>o remains confined in (— oo,mi) up to time t; we have 
Tl = for all i e {1, . . . , n} and then 1 — e"^'^''^'^ = 0. So, we can write that 

E^l - e-<'^^T,>^^^ ^ ^yy^y = e-<>'''^^>,Xt € dy,T,, < t)/dy 

= (l - e-<'^^^'--i > , Xt-r^^^ o e dy, < t) /dy 

= / E,{r„, eds}E„,(l-e-<'''T'->,Xt_, edy)/dy. 
Jo 
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Taking the Laplace transform with respect to i, this yields 



nOO 

/ e-^* [E^l -e-<'^^T'>,Xt e dy)/dy] dt 

JQ 



,-At 



<fJ..Tt> 



e r^{T„i e dt} /" 
Jo 

= (c-^^"i ) , y) - iP^ Ju^ , y)] . 

This proves p.Sp and p.4p ensues by integrating p.5p with respect to j/ on M. ■ 

3.2 Differential equations 

We deduce from Eqs p.ip and p.2p the famous backward Kolmogorov equations. 
Proposition 3.3 The functions ip^ ^ and tp^ ^ solve the following differential equations: 



Xt e dy)/dy] dt 



i=l 
n 



i'xjx,y) - Sy{x). 



These equations read more explicitly as follows: for a; G R, 



^ I (A + V'a.^Is;) - 1 if a; G (uj,t;j) and « g 
^'^ \Xip^^{x)-l if a; e and j G {0, 1, . . . , n}, 



(3.6) 



and for x g R\{?/}, 

T^xi^^^^x.y) = 



(A + ^j) i^^ JyX, y) if a; e (mj, Vi) and i G {1, . . . , n}, 

A t/'a ^1(2^1 y) if a; G (w^, Mi+i) and z G {0, 1, . . . , ri}. 



(3.7) 



Additionally, we see by p.ip and the regularity properties of that (f^ ^ is difFerentiable on R. So, 
we get the following conditions at points Ui,Vi, i G {1, . . . , n}: 



(3.8) 



In the same way, we see by p.2p that the function x 1-^ ^{x,y) is differentiable on R\{?/}. So, we get 
the following conditions at points Ui,Vi, i G {1, . . . , n}: 



ii^xjuf,y) ^i;^Ju^ ,y), ^[Juf^y) ^ ^[Ju, ,y). 
Moreover, due to (|2.1I) and p. 21) . we have the following conditions at point y: 

''^-{y ,y) — of'^y ^y)^-'^iy)- 



(3.9) 



(3.10) 



dx 



Different cases for y must be distinguished: y G for a certain io G {1, . . . ,n} or ?/ G (vio,Wio+i) 

for a certain iq G {0, 1, . . . , n}. 

In the forthcoming sections, we solve Eq. p.6p with conditions p.8p and Eq. p.7p with condi- 
tions p.9p - p.l0p . Let us introduce a basis of fundamental solutions {a.i,bi} of the linear second-order 
differential equation 'Du(x) = (A + /ij) u{x) for any i G {1, . . . ,n} as well as a basis of fundamental solu- 
tions {c, d} of the equation 'Du{x) = Ait(a;). These functions are chosen such that a^, c are increasing and 
6i, d are decreasing. 
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4 Solving differential equations (13.61) and (13.71) 
4.1 Solving Eq. SSM with conditions SSM 

In view of the differential operator T>, we see that the form of the solution of p.6|) is 



aiai{x) + (3ibi{x) + for x € {ui,Vi) and i £ {1, . . . ,n}, 

A + fit 

7i c{x) + Si d{x) + — for x g {vi,Ui+i) and i E {0, 1, . . . ,n}. 
A 



We have to determine the unknown coefficients 70, (5o, /3i, 7i, . . . , a„, /3 Put Vi — , ,,'7 - ■ 

For large enough negative a; (so that x < ui), we have E^(e^-^^"i) = c(a;)/c(ui) and p.4p supplies 



, , /I / n\ c(x) 1 

^-» = (A-^->^))il7) + A 



("1) 

This clearly implies that Sq = 0. Similarly, considering ip^ ^{x) for large positive x, we see that 7„ = 0. 
Next the regularity conditions p.Sp at Uj and Wj yield 



(4.1) 



'ai ai{ui) + Pi bi{ui) ~ 7^-1 c{ui) - Si-i d{u.i) = v^, 
ai a[{ui) + Pi h[{ui) ~ 7i_i c'(ui) - 5i„i d'(Mi) = 0, 
tti ai(i;i) + Pi bi{vi) - 7i c(wi) - Si d{vi) = i^j, 

^ai a'^{vi) + Pi 5-(wi) - -fi c'{v.i) - Si d'(vi) = 0. 

Let us introduce the matrices 

^-(?:)' «'=0:)- ^'-(i)' "'(^' =(::!:! ^ig)' .^'S 

Notice that Mi{x) depends on A + /i^ and N{x) depends on A. The system (|4.ip can be rewritten into a 
matrix form as 

{M,{ui)A, - iV(u,)B»-i - i^^Co, 
\Mi(^;,)A, - 7V(^;,)Bi = i/iCo. 

We extract from the first equation of (|4.2p the relationship Ai — Mi{ui)^^ N {ui)Bi^i + ViMi{ui)^^CQ. 
Plugging this equality into the second equation of (|4.2I) . we derive the recursive identity for Bi 

= [N{v,y'M,{v,)M,iu,)-^Niui)]B,^i + u,[Nivi)-^ M,{v,)M^iu^)-' ~ iV(i-,)"']Co. (4.3) 

Setting 

= N(vi)-HUvi)M,{ui)-^N{u,), Q, = N{v,)-^Mdvi)M,{u,)-^ - N{v,)-\ (4.4) 
we write (|4.3p more concisely as 



B, = P,B,^i + i/,Q,Co. (4.5) 
By iterating this recursive equality, we get for i G {0, 1, . . . , n} 

B, = P,P,^i . . . PiBo + {U,Q^ + U,^lP,Q^-l + l^^-2P^P^-lQ^-2 + ' ' ' + ^^lP^P^-l ■ ■ ■ P2Ql)Co 

with the conventions that PiPi^i . . . P2 ^ O and PiPi-i ... Pi = / if i = 0, and PiPi^i . . . P2 = I \i i = 1- 
The condition i^o = gives Bq = 70C0. Put, for i € {0, 1, ... , n}, 

R, = P,P,_i ... Pi, S,^ U,Q, + J/,-lP»Q,-l + V^-2P^P^-lQ^-2 + ' ' ' + • • • (4.6) 

By the aforementioned conventions, we have Pq = I, Sq = O, Si = i^iQi. With these settings at hand, 
we get an simple representation for P^: 

= [7oi?. + 5.] Co. 
Now, we use the condition 7„ = 0. By observing that 

7„ - (1 0)P„ = 7o(l 0)P„Co + (1 0)SnCo, 
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we deduce the value of the coefficient 70: 



70 



(1 0)5„Co 



(1 0)RnCo' 

Let us point out that (l 0)i?„Co (resp. (l 0)S'„Co) is nothing but the first entry of the matrix i?„ (resp. 

Sn). 

Finally, we express Ai by means of Bi: due to (|4.2|) . we have 

A, = M,{vi)-^N{v,)B, + u,M,{vi)-^CQ 

= [^oM,{v,y^N(vi)R, + M,{v,r^N{vi)S, + v,M,{v,)-^] Co- 
We sum up the results obtained in this section in the statement below. 

Theorem 4.1 The iterated Laplace transform of the probability distribution ofTt is given by 

1 



,{e-<''''^^>)dt = 



(1 0)Ah{x)A, 



A + fii 



with, for i G {1, . . . , n}, 

A,=M,{v,)-'N{v,) 
forie {0,1,..., n}, 



1 0)N{x)B, + - 
A 

(1 o)sJl 



for X e [ui, Vi] and i £ {1, . . . ,n\ , 
for X e {vi, Ui+i) and i G {0, 1, . . . , n\ , 



5. 



Ri 



Bi — Si 



1 0)i?: 

(1 Q)Rn(l 



A(A 



or 



imply 



where the matrices Ri, Si, i G {0, 1, . • . , n}, are defined by ^4-4^ |^.6p . 
Remark 4.2 The relationships i?„ = [joRn + Sn]Ca and _B„ = ^n^'o with Dq = 

S„ = (0 l)Bn = 70 (0 l)RnCo + (0 l)5„Co 
which gives the following expression of the coefficient (5„: 

<5„ - (0 l)SnCo - "1^"^" (0 l)i?„Co. 

Another expression can be derived by reversing the sense of the algorithm we used for solving the sys- 
tem Indeed, let us rewrite (H3t as B.^^i = P^^B,, ~ u^P~^Q^Co. Set 

P^ = Pr^ = N{ui)-^M,{u,)M,{vi)-^N{vi), Q, - -P^^Q^ - N{ui)-H'h{u,)M,{vi)-^ - N{ui)-\ 

The matrices Pi and Qi can be easily deduced from Pi and Qi by interchanging Ui and Vi. We have the 
following recurrence: 

which can be successively iterated from i = 1 to n; this yields 

-So — RnBn + •S'nC'o 

where 

i?„ = A A . . . Pn, Sn = VlQl + V2P1Q2 + VZP1P2QZ + • 

Since = 7oCo and i3„ = SnDg, we get 

7oCo = SnRnDo + S'nCo 

and we extract the simple expression 

^ (0 l)5„Co 
(0 l)i?„i?o' 



^nPlPl ■ ■ ■ Pn-lQn- 
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4.2 Solving Eq. ( l37fl) with conditions ( ISTOl ) and ( ISAOl) 

We have several cases to consider depending on the location of y inside or outside the intervals of the 
set E. 

First case: y £ [ui^ , Vig] for a certain io G {1, . . . , n} 
The form of the solution of (I3.7P is 



' ai{y) ai{x) + Piiy) h{x) for x G (w^, w^) and i € {1, . . . , n}\{io}, 

7i(2/) c(a;) + (5i(?/) ci(x) for x e {vi, Ui+i) and i G {0, . . . , n}, 

(y) ) + Ao,i(y)^Jo(a;) for x e {uto,y), 

(y) ) +/3io,2(y)&»o(a;) for a; e 



If y = (resp. Vi^), we consider that the interval {ui^, y) (resp. (y, Vi^)) is empty. We have to determine 
the unknown coefhcients 7o(y), (5o(y), (y), 7i(y), (^i (y), ■ • • , a»oa(y)' /3io,i(2/)> aio,2(y), Ao:2(y), ■ • • , 

an(y),/3™(y),7n(y),<5«(y)- 

For large enough negative x (so that x < min(j/, ui)), we have p^{x, y) — c{x) d{y) and E^; (e~^'^"i ) = 
c{x)/c{ui). Thus (I3.5P gives 



c{x). 



This clearly implies that (5o(y) = 0. Similarly, considering ij:^ {x,y) for large positive we see that 

7n(y) = 0. 

Next, conditions p.9p at Ui,Vi, i G {1, . . . ,n}\{io}, yield 

' ai{y) ai{ui) + fii{y) h{ui) ^ li^iiy) c{ui) + (5i_i(y) d{ur), 
ai{y) a',{ut) + My) h'^{ui) = li-i{y) c'{ui) +S,^i{y)d'{ui), 
ai{y) ai{vi) + ^y) bi{vi) = ji{y) c{vi) + Si{y) d{vi), 
a^iy) a'^M + A(y) K{v,) = 7,(y) c'(w,) + S,{y) d'{v,). 

Similarly, conditions p.9p at , give 

'aioi(y) ) + Aoi(y)^»o(w«o) = 7»o-i(y) c{uig) + 5ig^i{y) d{uig), 

"»o2(y) ) +/3^o2(y)^^o(^^io) = 7*0 (y) c( Vig j + Gig 

(y) d{v^g), 

.aig2{y) a'ig(vig) + I3ig2{y) Kg{vio) = 7»o(y)c'( Vig ) + Oig 
Additionally, conditions p.lOl) at y yield 

fa*oa(y) aio(y) + ftoa(y) ^«o(y) - aio,2(y) atoiy) - MAy) Kiy) = o, 
[atoAv) <(y) + fto,i(y) ^jo(y) -"io,2(y)<(y) - Ao,2(y)&jo(y) = '^(y)- 

Using the matrices introduced in the foregoing subsection and setting also 



we rewrite all these equations into a matrix form: for i G {1, . . . , n}\{io}, 

(M,{u^)A{y) = N{u,)B,^i{y), 
{M,{vM^iy) ^ Niv,)B,iy), 

for i — iQ, 

{M,g{U,g)A,gi{y) ^ iV ( W , „ ) „ _ J ( ^ ) , 

\M,g{V,g)A,g2{y) ^ N{V,g)B,g{y), 



(4.7) 



(4.8) 
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and finally 

M,„{y)A,,2{y) - M,,{y)A,,,{y) = -K{y)Do. (4.9) 

We extract from (jiT)) the relationships -Bj(y) = Niviy''- Mi{vt)Ai{y) and Ai{y) = Mi{uiy'-N{ui)Bt^i{y) 
which entail for i € {1, ... , n}\{io} 

B,{y) = P'A-i{y). (4.10) 
Therefore, since 5o{y) = and then i3o(y) = 7o(2/)C'o, 



B^{y) 



with, in view of (14.81) and (|49 



7o(2/)P,F.-i . . . Pi Co if i e {0, 1, . . . , io - 1}, 
PiFi_i . . . Pig+iBig{y) if i e {io + 1, . . . 



^..(y) = A^(«.o)"'^^.oK)^2o2(2/) 

= N{v,„)-^M,M[A^ol{y) - «(2/)Af,„(y)-ii?o] 

= Niv,„)-Hl,,{v,,)[M,„{uJ-'N{u,,)B,,^i{y) - ^,iy)M,,{y)-' Do] 

= 7V(t;,J-iM,„(«,J[7o(2/)M,„(u,J-i7V(M,jF,„_iP,„_2 . . .PiCq - k(2/)M,„ (y)-ii^o] 

= 7o(2/)P,:„P»o-i ■ • --PiCo - «(y)iV(,;,J-iA4,(i;,jAf,„(y)-iZ?o. 

As a byproduct, recalling that we set Ri — PiPi^i ... Pi and putting, for i e {0, 1, . . . , io — 1}, 

U, = P,P,-i . . . P,„+iN{v,„)-Hd,M (4.11) 

with the conventions that PiPi-i . . . Pig+i = I ii i ^ iq and PiPi^i . . . Pi^+i = O if z < zq, we obtain the 
following expression for Bi{y): for any i g {0,l,...,n}, 

B^iy) = loiy)R^CQ - Kiy)U,M,^{y)-'Do. (4.12) 

Now it remains to determine "fa{y). For this, we invoke the condition "fniy) = 0. Since 

7„(y) = (l 0)P„(y) =7o(2;)(l 0)P„Co - «;(2/)(l 0)C/„M,„(y)-iDo, 

we deduce the value of the coefficient 70(2/): 

^ (1 0)U,,M,,{y)-'Do 

{T^l^o — • 

Finally, due to ()4.7p and ()4.12p . we can express Ai{y) by means of Bi{y), for i e {1, . . . ,n}\{io}, as 

A{y) = Mdvi)-^N{vi)B,{y) = Kh{vi)-^ N {v,) [^o{y)R,Co - K{y)U,M,,iy)-' Do]. 
Moreover, by (g^), (gS]) and (j4l^ . 

^2ai(y) = A^2o(w2o)"^^K0S2o-i(y) = 7o(y)M,„(u..o)"^^K)^.o-iCo, 

^2o2(y) = M,„{v,„y^N{vM,{y) = 7o(y)M,„(u,J-iiV(u,jP,„_iCo - K{y)M,,{y)-^Do. 



Remark 4.3 As in Remark 14.21 we can obtain two expressions of the coefficient (5„(y). Because of the 
relationships Bn{y) = jo{y)RnCo - K{y)UnMig{y)^^ Do and P„(j/) = Sn{y)Do, we can see that 

J„(y)^.(y) "1^"^^"^^^''^° (0 l)P„Co-.(y)(0 l)t/„Af,„(y)-ii^o. 

Another expression can be derived by reversing the sense of the algorithm we used for solving the sys- 
tem (|i?7| . Indeed, let us rewrite (03]) as Pi-i(y) = PiBi{y) for i e {1, . . . , n}\{io} (recaU that P,; = P^^). 
We have 

Po(2/) = Pi . . .P,„_iP,„_i(2/), P,„(y) = 4+1 . . .P„P„(y). 
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Moreover, by ()4.8p and (|4.9p . we successively have 

- 7V(7.,„)-iAf,„K)[M,„K)-i7VK)B,„(y) + «(y)Af,,(y)-ii5o] 

= hMv) + K{y)Niu,,y'M,„{u,„)M,„{yy'Do 

= Ao/'.o+i • ■■PnBniy) + K{y)N{u,„)-'M,^iu,,)A'l,^{y)-'Do 

and then 

Bo{y) = A . ..PnBniy) + n{y)Pi . . . h„^iN[u,„)-H'h,{u^,)M,,{y)-^ D^. 
We rewrite this last equality as 

7o(y)Co = Sniy)R„Do + niyWoM^M'^ Do 
with Rn — Pi . . . Pn and JJq — Pi . . . Pig^iN (uig)^^ A4ig{uig) . We finally get the following expression of 

(0 l)UoM,^iy)-^Do 



(0 l)RnDo 



Second case: y g (vig , Wio+i) for a certain iq G {0, . . . , n} 

Since the computations are very analogous, we briefly outline the corresponding matrix equations. 

Set 

The equations write as follows: 

M,{u,)A,{y) = N{u^)B,^i{y) for i e {1, . . . ,n}\{io + 1}, 

Mi{vi)Ai{y) = A^(v,)B,(2;) for z e {1, . . . , n}\{io}, 
M^MA„{y) = N{v,g)B,,iiy), 
Mig+i{u^g+i)A^„+i{y) = N{u^„+l)B,g2{y), 
[Niy)B,M - N{y)B,M = -K{y)Do. 

We deduce that 

{io{y)P.P.-i ■ ■ ■ Pi Co if z e {0, 1, . . . , io - 1}, 

\PiP,^i . . . P,^+2Bta+i{y) if i e {io + 1, ■ ■ ■ 



B^iy) 
with 



Bio+i{y) = iV(ui„+i)"iM,„+i(t;,,+i)A,„+i(y) 

= iV(t>i„+i)-iM,„+i(i;,„+i)Mi„+i(ui„+i)~iiV(ui„+i)Bio2(y) 
= Pio+iBio2{y)- 

We need to compute Big2(jj)'- we have 

i?.o2(y) = S.oi(2/)-«(2/)^(2/)"'^o 

with 

i3,„i(y) = N{v,,)-'M.,,{v,,)A,„iy) = N{v,,r^Nh,{v,,)M,,{u,,)'^N{uMo-i{y) 
= P,,B,,_M = 7o(y)P,oP.o-i ■ ■ • PiCo. 

Thus 

= [7o(tj)P.o • • • ACo - n{y)N{y)-^Do\ 

and for i G {ig + 1, . . . , n}, 

= 7o(2/)P»i'^-i • • • -Pi Co - K{y)P,P,_i . . . P,„+iN{y)-'Do. 
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Setting, for i G {0,1,..., n}, 

=P,P,_i...P,„+i, (4.13) 

with the conventions that Vi — I ii i = iq and Vi = O ii i < io, we obtain the foUowing expression for 
Bi{y): for any i G {0, 1, . . . , n}, 

B^{y) = My)R^CQ - K{y)V,N{y)-'Do. 
The condition ^n{y) = supphes the value of the coefficient 70 (y) within Bi{y): 

(1 Q)V^N{y)~^Do 



(1 0)i?„Co 



Remark 4.4 As in Remark 14.31 we can obtain two expressions of the coefficient 5n{y)- We only outline 
the main details. On one hand, the first expression writes 

5n{y)^<y) ?/"fi^l,'^° (0 l)i?„Co-«(y)(0 l)Vr.N{y)-'D^. 

On the other hand, we have 

B^{y) = Pi . . . Ao-iS,:„-i(y) = Pi . . . P^oB^oliy) = A . . . Ao[S.„2(y) + «:(2/)Af(2/)-^i?o] 

= Pi . . . Ao+iS,„+i(y) + k(j/)Pi . . . P,„N{y)-^Do = A • • ■ P„P„(y) + k(2/)Pi • • ■ P^,N{y)-^DQ 
which we rewrite as 

lo{y)Co = Sniy)RnDo + n{y)VoN{yr^Do 
with Vq = A • • • Pjo- We get the second expression of Sn{y): 

(0 l)VoNiy)-^Do 



1 P„i?o 



Finally, we can express Ai{y) by means of Bi{y): for i G {1, . . . ,?T.}\{jo}i 

= M,{vi)-^N{v,)B,{y) = M,(«0-^Af(vO [7o(y)i?.Co - K{y)V,N{y)-^ D^] 

and, for i = io, 

A,„{y) = M,„K)-i7VK,)P,„_i(y) = A/,„ K)-17V(m,J [7o(y)P.„_iCo - K{y)V^,.iN{y)-^ Do]. 
We sum up the results obtained in this section in the following statement. 

Theorem 4.5 The iterated Laplace transform of the joint probability distribution of (Tj, Xt) is given by 
the formulas below. 

1) If y & [mjo ! ■^io] /o'" o, certain io £ {1, . . . ,n}, 

(1 0)Mt{x)At{y) for X e [ui,Vi] and i e {I, . . . ,n}\{io}, 
(1 Q)N{x)Bi{y) for x G (t;^, Ui+i) anrf i G {0, . . . , n}, 



e-^* [E4e-<'^^T'>,Xt G dy) /dy] dt 



(1 0)M,„{x)A^„i{y) for X e [u,g,y], 
(1 0)M,o(x)Ai(,2(y) /or a; G [y,t;io]; 



with, fori G {1, . . . , n}\{io}, 



(1 o)c/„Af.o(y)-'(J) 



1 P„ 



r,.{q]-u,m,m 



(1 o)c/„Af.o(y)"'(5 



^»o2(y) = K(y) 



(1 0)P, 



(1 0)C/„Af,„(y)-i('j 



(1 0)P„ 
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and, for i € {0, 1, . . . , n} , 



(1 0)UnAhM-'(^^ 

Rii „ 



(1 0)i?„( J 



1^0 



where the matrices Ri,Ui are defined by ^-4^ , |4.6p artrf j 
2)Ifye (wio, Uio+i) /or a certain io e {0, . . . , n}, 



(l 0)A/i(a;)Ai(y) /or a; e [u^, u^] anrf i G {1, . . . , n}, 
_ J (l 0)N{x)Bi{y) for X e {vi,Ui+i) and i £ {0,. . . ,n}\{io}, 



1 0)A^(a;)Bi„2(y) /or a; e [y, 
with, fori G {1, . . . , n}\{io}, 



(1 0)KA^(,)-(J) 



(1 0)F„iV(,)-(j) 



and 



(1 0)KA^(2;)-i(j) 



(1 0)i?„ 



B^o2{y) = K{y) 



(1 0)KA^(2/)-i(; 



(1 0)i?„(^J 
(1 0)KiV(2/)-i(^j) 

(1 0)Rn(l 



where the matrices Ri,Vi are defined by J^.-^p , |^.6'| j anrf J^. 

5 Examples 

5.1 Case of one bounded interval 

In this part, we focus on the set E made of one two-sided interval E — [u,v]. This classical case 
corresponds to the values of the parameters n ^ 1 and ui = u, vi = v, = Tt. We relabel the argument 
fii into /i, the functions ai, hi into a, b and the related matrix Mi into M. The settings write here Rq = /, 
5*0 = O, i?i = Pi = N{v)-^M{v)M{u)'^N{u), Qi = N{v)-^ M{v)M{u)-^ - N{v)-\ Si = j^Qu 

i?i = A = N{u)-^M{u)M{vy^N{v), Qi = N{u)-^M{u)M{v)-^ - N{u)-\ Si = Qi. 



Probability distribution of Tt 
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The distribution of Tj is characterized by 



Jo 



(l 0)N{x)Bo + - for X e (-00, u], 

A 



(1 0)M{x)Ai + 
(1 0)N{x)Bi + 



A 



for X e [u, t;] , 
for X e [u, +00). 



Observing that M{vy^N{v)Pi = M{u)-^N{u) and M{v)-^N{v)Qi + M(i;)-i = M{u)-\ the matrix 
j4i can be simplified into 



A,=,oM(.)-iV(.)(j)+3^MH-(;) 



and the other matrices are given by = 7o fg j and Bi = 5^ f M with 



(1 0)Qi 



70 = - 



'^'^'^ (1 o)pJl 



Si = - 



(0 i)Oi 



(0 i)a(; 



Moreover, 



(1 0)iV(a;)Bo = 7oc(a;), (l 0)N{x)Bi = 6id{x). 



Joint probability distribution of {Ti, Xi) 

1) Assume that y € (— cxd,u]. This case corresponds to = and we have Vq = I, Vi = Pi- The 
distribution of (T(,Xt) is characterized by 



POO 

/ e-^* [E, (e-'^^* , Xt e Ay) /dy] dt = < 
Jo 



{ (1 0)iV(x)Boi(y) for (-00,2/], 

(1 0)iV(x)Bo2(y) fora;e 

(1 0)M(a;)^i(y) for x e [u,w], 

. (1 Q)N{x)Bi{y) for a; e [v, +00). 



Since M(t;)-i7V(u)Pi = M(M)-iiV(M), we have for 

A^{y)=^o{y)M{v)-^N{v)(^^-K{y)M{u)-'N{u)N{y)-^(^^, 
and the other matrices are given by 

Soi(y) = 7o(2/)Q, So2(y)=7o(y)Q 

-1^0 



Biiy) = Pi 



My)[Qj -K{y)N{y) 



with 



(1 0)PiNiy) 



70 (y) = K{y) 



(1 o)P.L 



2) Assume that y € [uj-y]- This case corresponds to io = 1 and we have Uq = O, Ui = N{v) ^M{v). 
The distribution of {Tt,Xt) is characterized by 

{{1 0)N{x)Bo{y) for a; e (-oo,tt], 
(lO)M(.Mu(,) for. e 
(1 0)M{x)AM forxe [y,^;], 
(1 0)N{x)Bi{y) for X G [v, +00). 
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The matrices are given by 

Aiiiy) = ^o{y)M{u)-'N{u) Q , AM = joiy)M{u)-'N{u) Q - «(y)M(y)-i (^f^ , 

Boiv) = 7o(y)Q, = 7o(2/)PiQ - «(y)t/iM(y)-i(^J 

with 

(1 0)[/iM(y)'i[j 

70 (y) = K{y) 



(1 0)P, 



3) Assume that y € [v,+oo). This case corresponds to io = 1 and we have Vq = O, Vi = I. The 
distribution of {Tt,Xt) is characterized by 



/ e-^' [E, (e-^^* , Xt e dy) /dy] dt={) 
Jo (I 



(1 (})N{x)Bo{y) for a; e {-oo,u], 

(l 0)M(a;)Ai(y) for a; e [u,?;], 

0)iV(.T)Sii(y) forxe 

(1 0)N{x)Bi2{y) forxe [2/,+oo). 



The matrices are given by 



My) = lo{y)M{v)-'N{v) (^gj , Bo{y) = 70(2/) (^gy , 
Biiiy) = 7o(y)Pi Q , = loiy)Pi Q - n{y)N{y)-' (^^ 

with 

(1 0)Niy)-4f 

70 (y) = 7- ^ 

(1 o)Pi(J 

5.2 Brownian motion 

In this example, we take for {Xt)t>Q linear Brownian motion rescaled such that E{X^) = 2t, that is 
Xt = V^Bt where {Bt)t>o is the standard Brownian motion satisfying K{Bf) = t. This choice which 
corresponds to a{x) = \/2 and t{x) = is done for simphfying the forthcoming settings. We take for 
a, 6, c, d the functions 

a{x) = e^'^^, b{x) = e-^^"", c{x) = e^^, dix) = e"^^. 

The potential writes p^{x,y) = ^^e~^l^~^l and then w{x) = 2-\/A and k{x) = 1. The related 
matrices M and A'' write 



M{x) = 

N{X) : 

and their inverse are given by 

M{x)-^ 

N{x)-^ 



1 fy/XTJie-^^' 



2VATm V \/ATMe^^+''^ _eVA+^x 
1 /v^e"^^ e-^^\ 



2\/A V \/Ae^^ -e^^ 
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We have 



M{y)M{x)-^ 



sh(v^AT7i(y - x)) 



sinh(vX+7t {y - x)) ' 

y^/X + ^ sinh(VA + fx{y - x)) cosh(i/ X + ii{y - x))^ 
sinh(^)' 

cosn(z) 



yx+T^ I, 

^y/X + fj, sinh(^) cosh(2) 
where we set z = y/X + jj. {y — x) for Ughtening the text. Next, 
N{y)-^M{y)M{x)-^ 



cosh(z) + ^/^sinh(z^ 



cosr 



sh{z)-^smh{z) 
Making use of the elementary identity 



cosh(z) + JjT- sinh(z) 



cosh(2:) + sinh(^;) 



cosh X + a sinh x = 1 + 2 sinh^ — + 2a cosh ^ sinh ^ = 1 + 2 sinh ^ cosh — + sinh ^ 



we get 

7V(y)-iM(y)M(a;)-i - N{y)-^ 



-Vxv 



sinh(f) 



y^cosh(|) + sinh(|) 



^cosh(|)+sinh(|) 



-VXy 



y3^cosh(|) + sinh(i 



^ cosh(f) - sinh(f ) J 



Now, 

Ar(y)-^M(j/)M(a;)-^Af(a;) 



/ pVX (x-y) 



2 COS 



X 
A+/i 



A 



^ sinh(z) 

e^^(^-) [2cosh(^) - yjl) sinh(^) 

Observing that 

a H — I sinh x = 2 ( cosh^ — + sinh^ tt ) + 2 ( a H — ) cosh — sinh — 
a J \ 2 2J\aJ 22 

2 / , X . , a;\/ , a; . , a;\ 
= - ( a cosh — + smh — 11 cosh + a smh — 1 , 

Qi \ ^ Zi / \ Zi Zi J 

we obtain 

N{i))-^M(%))M[xY^'N[x) 

1 /e^^^-s') [\/Acosh(f ) + \/ATMsinh(f )][VATMCosh(f ) + \/Asinh(|)] 
v/A(A + m) \^ -/ie^(^+f)cosh(f)sinh(f) 

/xe-^(^+2') cosh(|)sinh(|) 
ev^(!/-x) [^/Acosh(l) - ^/ATMsinh(f)][^/ATMCosh(f) - \/Asinh(f)]^ 
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5.3 Case of one bounded interval for Brownian motion 

We now consider the sojourn time of Brownian motion in the interval [u, ?;]. In order to lighten the 
we only compute the distribution of Tj, that of {Tt,Xt) being more cumbersome. 
For evaluating 70, we compute, with w = |\/A + /j, {v — u), 

(l\ g\/A(«-t)) 
= [vAcosh(w) + -v/A + /Usinh(«;)][VA + /iCOsh(«;) + VXsmh{w)\ 
x/A(A + u) 



and 



VA(A + /x) 

(1 0)(5i y^j = sinh(u') [^/A + /xcosh(w;) + vAsinh(tz;)] 

from which we deduce 

^ (1 O)Qi(o) ^ e-^"sinhH 

70 = -■ 



A(A + /i) o)Pi^'^^ A^XTm VAcosh(w;) + VAT7tsinh(w;) " 
For evaluating ^i, we compute 

(q\ gVA(M— f) 

) = —j^=^= [\/Acosh(u;) + ^/\ + n sinh(«;)] [-\/A + cosh(«;) + \/Asinh(w)] 
^/ a/A(A + u) 



and 



= — ^ sinh(M;) [^J \ + 11 cosh.{w) + VAsinh(M;)] 



from which we deduce 

(0 l)Q,{\ 



di = -■ 



^ e^''sinh(w) 



X{X + lj) I'^p^^^ AVA + n ^/Xcosh{w) + + iismh{w) ' 
This yields 

'l\ _ ^ sinh(w)e^(^-") 



(1 0)Af(a;)Bo=7o(l 0)Ar( 



Ay'A + /X •\/Acosh(w) + ^/A + ^ sinh(w) 



(1 ,)Ni.)B,=5,{l 0)iV(x)r?U--^ ^ sinhHe^-) 

^ ^ ^ ^ Vv AVAH^ VAcoshH + VATTIsinhH 

On the other hand, 



MW-(j)=ie 

M(u)-^7V(u) 



1\ _ c^" /^XTTIe-^^" e-^^"Y 1 
Oy ^ 2VATm V yXTTIe^'^" _eV^"Ji^VA 

e^" /[yXTTt + yA] e-^^^"\ 
~ 2^/ATm V [\/aTm - \/A] e^^" y 
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and then 



^1 =7oM(w)-i7V(u)Q + 
/ 



-M{u) 



2X{X + fi) 



1 - 



A(A + /i) 
(VA + /i + VA) sinh(u;) 



\/Acosh(u>) + ^/X~\~J^ smh{'w) 
+ fi — a/a) sinh(t(;) 

V L ■\/Acosh(w) + a/A + ^ sinli(w) 



1 - 



[cosh(w) - sinh(w)] e~^-^+^" 
2\/A (A + ^) \/A cosh(w) + VA + ^sinh(w) V [cosh(w) + sinh(u;)] e^/^" 

-VA+7r(«+ti)/2 



2VA (A + fi) VAcosh(w) + VAT7Isinh(w) V cv^("+")/2 



Now, 

(1 0)M(x)Ai = 



/c-\A+m(«+^')/2\ 



2\/A (A + /i) \/Acosh(w) + -y/A + fismh{w) 

/i cosh(VA + ^ (a; - (m + w)/2)) 
\/A(A + ^) \/Acosh(u') + ^/A + /isinh(w) 

We sum up below the results we have obtained. 

Theorem 5.1 The distribution of the sojourn time Tt inside the interval [w, v] for Brownian motion is 
characterized by 



f 

Jo 



e^^*E^(e~^'Mdi = <! 



'1 


1 - 
1 


A 


A 


-f /i 


1 


1 - 


^ A 



/i sinh(7«)c^(^-") 



\/X + j-i v'Acosh(u>) + yO^T+Ji smh{w) 

^ l_i cosh(\/A + fi{x - {u + v)/2)) 

\/A \/Acosh(?i;) + v'A + /Lisinh(u;) 

fi sinh(w)e^(''-=^) 



^/X + fj. -\/Acosh(w) + y/X~r]lsmh{w) 



for X G {—oo, u], 
for X e [u, v], 
for X € [v, +oo), 



where w = ^\/A + /i (w — u). 

We retrieve the well-known distribution (1.7.1), p. 140 of [I]. 



5.4 Local time in a finite set for Brownian motion 

Let us apply our results to the following set (ui, . . . , u„ are real numbers such that ui < • • • < m„): 

n 

= IJK - e,Ui + e] 

i=l 

where e > is subject to tend to and denote Tt ^ = 1e- {Xg) ds. Set also uq = —oo and = +oo. 
The local time in the set {ui, . . . , u„} of Brownian motion {Xt)t>a up to time t is defined by 

Lt = lim - Tt.e- 

e-i-O e 

Of course, as previously, we can decompose Lt into the sum 

n 
i=l 

where, for any i g {1, . . . , n}, 

1 /■* 

ij = lini-/ l[„._j.„^+j](Xs) ds. 

e->-0 e Jq 
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Let us introduce the vector of local times at each point u;: Lf = (L], . . . ,L"). Set, for i €E {1, . . . , n}, 

P,,e = N{u, + e)-^M,{ui + e)M,{ui - e)-^N{u, - e), 
Qi^e = N(u, + e)-^hh{ui + e)M,{ui - e)-^ - N{ui + e)-\ 

Ri.e = Pi,ePi-l,e ■ ■ ■ Pl.s, 

Si.e = l^i,eQi,e + ^i~l,ePi,eQi~l,e + t^i~2,ePi,ePi~l.,eQi-2,e + ' ' ' + Vl,ePi,ePi-l,e ■ ■ ■ P2,eQl,s, 
J^i,s — ^i,e'~^Q — J- — TTj: —J^i,e'-^Q- 

As usual, we set i?o,e — I, = O and So.e = — o)fl"'^c° C'o- 

In this part, we compute the following limit, for x G [ui, u^+i) (which entails that x € (ui +e, w^+i — e) 
for small enough e) and i S {0, 1, . . . , Ji}, 

/ e"^* E^(e-<'^'^*>) At = lim (l 0)iV(a:)B,,e + - 



where 



7V(a;) 



e" e 
%/Ae^^ -%/Ae-^^, 



Set e = e^/\ + ]jLi/e. Using the results of Section [5?^ we have for j g {1, . . . , n} 

I / e-2ev^[v/Acosh(e) + + /Aj/esinh(e)][^A + ^i/ecosh(e) + v/Asinh(e)] 

VK>^^nMje) y -(^,/e)e2"'^cosh(e)sinh(e) 



(^^/e)e-2"'^ cosh(e) sinh(e) 
[•\/Acosh(e) — a/A + /ii/£sinh(e)][-\/A + Hi / e cosh{e) — \/Asinh(e)]^ 



and 



Qi.e = sinh(e) 



(tti+s) 



A 



cosh(e) + sinh(e) 



cosh(e) + sinh(e) 



„v/A(,.i+e) 



1+77-71 cosh(e) + sinh(e) 
cosh(e) - sinh(e) 



By the elementary asymptotics for e — > 0+ 

VAcosh(e) + \/\ + ^i/e sinh(e) ^ \/A + ^i, 
a/a + /Ui/e cosh(e) + \/Asinh(e) ^ \J [iije, 
cosh(e) sinh(e) ^ ^/Jki, 

we get, for i e {1, . . . , n}, lim5^o+ Pr,e = A, lim£^o+ Qi,e = Q%, linij^o+ i?i,e = Ri, lim£_).o+ S'i,e = 5^ 
with 



1 / VA + /x^e' 
VAl-A^.e2^'" TA- 



-2\/A«i 



Va 1^- 



0' 



Ri — PiPi-l ■ ■ ■ Pi, Si — — [Qi + PiQi-l + PiPi-lQi-2 + • • • + PiPi-l ■ . ■ P2Q1] • 



We also have, for i e {0, 1, . . . , n}, limj^o+ Pi,e = Bi where 



Bi = Sil ^] - joRi ( Q I with 70 



(1 0)Sn 



(1 0)Rn 



(5.1) 
(5.2) 

(5.3) 



We can now state the following result. 
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Theorem 5.2 The iterated Laplace transform of the vector of Brownian local times lit o,t points mi , . . . , u„ 
is given, for x £ (u^, iti+i) and i G {0, 1, . . . , n}, by 



/ e-^* E,{e-<>''^'>) dt (1 0)N{x)B, + - 
Jo ^ 



A 

where the matrices Bi, i G {0,1,..., n}, are given by i5.j)) . i f 5. ^|) and kS.S^) . Additionally, this formula 
holds also for x — Ui, i G {l,...,n}. In particular, the iterated Laplace transform of the Brownian 
local time Lt in {ui, . . . , Un}, namely e^^* E^; (e^^^*) dt, can be deduced from the previous formula by 
choosing fi — (fi, . . . , /i). 

Proof 

It remains to prove the assertion concerning the case x = Ui. Observing that we have 



and 



(1 0)iV(..)Q. = ^ (e^- e-v^-)(!.'eSr. o) 



(0 0), 



we deduce that 

(1 0)iV(M,)i?, = (1 0)N{u,)P,P,.i...Pi^ {1 0)iV(M,)A-i...A = (1 0)N{u,)R^-i 

and 

(1 0)N{u,)S,^j[{l 0)N{u,)Q, + {l 0)7V(MOA(g^-i + ---+A-i...AOi)] 
= i(l 0)7V(mO(Qz-i + --- + A-i...P2Qi) = (1 0)N{u,)S,.i. 

As a result, since Bi is linear combination of the matrices Ri and S'i, we have 

(1 0)Niu,)B, = (1 0)N{ui)B,-i. 

This proves that 

hm / e-^*E^(e-<''^^*>)dt= lim / e~^* E^, (e-<''^^*>) dt 

x^uf Jo x^u^ Jo 

and then, by continuity with respect to x, 

poo pOO -i 

/ e-^*E„^(e-<'*^^'>)di = hm / e-^' E^{e-<''''^'>) dt ^ {l 0)N{u,)B, + -. 

■ 

We end up this part be considering the particular cases n — \ and n = 2. 
Local time in 

Suppose that n — 1 and set ui = u, fii = fi. We have Ro = I, So = 0,Ri = Pi,Si = jQi- Therefore, 
^i-yxl-Me^^" VA-/J' "^1- A3/2 l-cv^" oi' 



Ho 



/ie-^" /I 
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and 



Next, 



(1 0)5 



(1 0)i?,^Q 
1 0)N{x)Bo 



1\ '\0J A3/2 



a(Va + m) 



(1 0)N{x)Bi 



g%/A(«-an) 



A(VA + /i) A(VA + /i) 

As a result, the iterated Laplace transform of the local time Lt in {u} is given, for any x S M, by 

/ e-^*E4e-^^')d<= i- l-^=^^c-^l^-"l . 

We retrieve formula (1.3.1), p. 126 of [T]. 
Local time in {u, v} 

Suppose that n — 2 and set ui — u,U2 — v, /ii — fi, ^j,2 — v. We have Rq — I,Sq — 0,Ri — Pi, S 
iQi,i?2 = P2A,52 = i(<32 + AQi),7o= [(1 0)^2 Q] /[(I 0)i?2Q]. Explicitly, 



0' 



' yXl-Me^^" VA-/ii' ' A3/2 Ue^" oi' 



-R2 = ^ 



1 / %/A + z/ i/e"2v^"\/ VA + Ai /ie-2\^"\ 
A l^-t/e^^" VA-i/ yl^-^e^^'' VX- fi ) 

if (VA + ^)(VA + i/) -Mi^c2^("-^) (VA + t^)MC"2^" + (\/A-^)i/e-2^^' 
A i-(yA + z/)Aie2v^«- (VA + //)z.e2^'' (VA ~ - f^) - M^^e^v^^"""^ 



Hence 



( e^^^ 










A2 \ 



2-2v^A/e-^" 0^ 
fX-v ji-e^" 0, 



J_ /j/(%/Ae-^^ -^e^("-2t))^ _^ (VA + iy)/ie-^" 0^ 
A2 I -i/(VAe^^-Aie^(2t--")) - (\/A-j/)/^e^" 0, 



70 



z^(VAe-^" -/xe^("-2^.)^ ^ (VA + j/)Aie-^" 
A[(A/A + Ai)(\/A + z/) -;uz^e2VA 

//[VA + 1/(1 - e^(«-"))]e-^" + i/[\/A + m(1 - e^("-^'))]e-^'' 
A[(VA + M)(yA + I/) - /ui/e2^(«-'')] 



Therefore, since = — 7o ( q ) , we get 



(1 0)iV(x)Bo 



[\/A + 1^(1 - e^("-''))]e^(^-") + i/[VA + m(1 - e^("-"))]e^(^-") 



Now, 



20 



Straightforward computations show that numerators of the entries of Bi can be simplified into: 



Thus 



and then 



Bi 



(VA + /i)(\/A + v)- /iiye2^("-'') 
(\/A + m)(\/A + v)- /iiye2^("-'') ' 



(1 0)Ar(a;)i?i 



'a[(\/A + m)(\/A + I/) -/iz/c2^ \^/i[VA + i^(l-e^("^"))]e^^ 

j/[VA + - e^("-''))]e^(^-'') + /Lt[VA + ;/(! - e^("-''))]e^("-^) 
A[(^/A + ^)(v^ + J/) - Aii^e2v^ 



Furthermore, 

_ ^ / zy(VAe-^'' -/le"^^""^"') + (\/A + i^)Aie-^" -7o[(\/A + /i)(\/A + iy) - ^z/e^^^ \ 
" A3/2 l^-i.(VAe^"-//e^(2.-«)^ _ (VA - j/)Aic^" + 7o[(\/A + i^)/ie2^- - (\/A + Ai)i^e2^"] J " 

Obviously, the first entry of B2 vanishes. Tedious computations show that the numerator of the second 
entry of B2 can be simplified into: 



Me^"+7o(VA-Ai) 



(\/A + yu)(\/A + v) - /ii^e2v^(«-") 



Thus B2 = -^2 ( ^ I with 



h = 



and then 



^1 0)N{x)B2 = - 



^[%/A + 1/(1 - e^(«-"))]e^'' + + fi{l- )]e^" 
A[(VA + /x)(VA + i/) -/xi/e2v^ ("-")] 

/i[V^+ i/(l - )]e^("-^) + i/[VA + /i(l - e^("-''))]e^("~"') 

A[(%/A + /i)(\/A + J/) - Aii^e2^ 



As a result, we can see that the iterated Laplace transform of the couple of local times (L", L^) at w 
and V can be expressed by the unified formula, for any a; e M, 



We retrieve formula (1.18.1), p. 150 of [1 . Consequently, the iterated Laplace transform of the local time 
L"'" = L" + Lj in {u, v} is given, for any a; e M, by 

[%/A + /i(l - [e-v^l^""l + e-VA 



-At 1 



.(e--r)di = i 
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